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1 Introduction 

Throughout, for readabihty, we often use shorthand vector notation such as t for 
{ti, . . . , td) and dt for dti ■ ■ ■ dt^- 

A d-parameter (standard) Brownian sheet B = (i?(t):t G [0,1]'^) is defined to 
be a real- valued Gaussian random field with continuous sample paths, mean func- 
tion Ei?(t) = 0, and covariance kernel Ko{s,t) := E[i?(s)S(t)] = W^-^^min{sj^tj} 
for s,t G [0, 1]'^. Notice that the covariance operator with kernel Kq is the tensor prod- 
uct of d copies of the covariance operator of Brownian motion. Let m = (mi , . . . , rad) be 
a fixed d- vector with nonnegative integer components. We define in-integrated Brownian 
sheet by 

(1.1) X^{t):= ... n ^'^ B{du,,...,dud). 

Jo Jo fJi W 

It follows immediately that Xm is a mean-zero Gaussian random field on [0, l]'^ with 
covariance kernel 

minis, ,t,} _ ^.^^mj _ ^.^m, 



and so the covariance operator of is simply the tensor product of the covariance 
operators with kernels , which are the covariance kernels of mj-integrated Brownian 
motions. Observe that Xq is simply d-parameter Brownian sheet B. 

Remark 1.1. To motivate the definition of an integrated Brownian sheet, recall that 
m-times integrated Brownian motion Xm can be defined (naturally enough) in terms of 
Brownian motion B by 

(1.2) X^{t):= / ••• / B{si)dsi...dsm 

Jo Jo Jo 

for integer m > 1, with Xq := B. It is not difficult to see that Xm has the same 
distribution as the process with value at time t given by 

(1.3) ^ [\t-urdB{u). 

ml Jo 

That is, the m integrations in H1.2() can be collapsed to the one in ()1.3|) . Similar 
motivation can be given for our definition of m-integrated Brownian sheet. 

Consider the squared L^-norm of X^: 

"1 rl 
r2 — T/2 / / '\r2 



(1.4) v^^y^,= ... xi,{t)dt. 

Jo Jo 

The classical Karhunen-Loeve expansion tells us that X^ has the same distribution 
as the process (X]n V^^^^n ' where the ^n's are i.i.d. standard normal random 
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variables and the ip^s form a complete orthonormal system of eigenvectors, with cor- 
responding eigenvalues On, for the covariance operator Am ■ ^^([0, 1]"^) ^ L^{[0,1]'^). 
The spectrum a{Am) = {a-n} is the product of spectra of the covariance operators 
of the associated mj-integrated Brownian motions: 

O-(^m) = 0-(-4mi) ® • • • O Cr{Ama), 

where represents elementwise set multiplication: S T := {st : s £ S and t G T}. It 
therefore follows that, in distribution, 

oo oo 

(1-5) = E • • • E «n^n, 

where On = an(m) = ani{mi) ■ ■ ■ an^{md) and {anj{mj)) are the eigenvalues for Amj- 

We are interested in deriving a "strong" small-deviations result, i.e., in computing 
the lead-order asymptotics as e ^ 0+ for the small-ball probability F{V'^ < e) (not just 
for its logarithm). In light of the representation H1.5() . this can be done using Sytaja's 
Tauberian Theorem |15j : 

Theorem 1.2 (Sytaja |15p. Suppose that an > for all n with ^ ^''^'^ 

that (^n)n>i CLi"e i.i.d. standard normal random variables. Then as e ^ 0+, 

P {fl^^^n <ej= (1 + o(l)) [-2n{x*fh"{x*)]~'^\xp{-[h{x*) - ex*]} 

where h(x) := — logE exp {— x ^^-j^ a„^^}, x > 0, denotes the log Laplace transform 
and x* = x*{e) is defined implicity in terms of e: 

(1.6) h'{x*)=e. □ 

In this paper we first show how to derive, for arbitrary m, a complete asymp- 
totic expansion as x — > oo for the log Laplace transform of V'^, namely, hm{x) := 
— logE exp{— xF^), and for each of its derivatives. Our method proceeds in two steps. 
First we study carefully the one-dimensional functions ■ Then we apply Mellin- 
transform techniques that prove to be quite powerful in handling our "tensored" pro- 
cesses. Although Mellin transforms enjoy common use in the analysis of algorithms |5, 
their application here is, to our knowledge, the first in studying small or large devia- 
tions. As the reader will see, their use is very natural and makes many computations 
completely transparent. 

We also show how one can develop the explicit asymptotic behavior of x* in terms 
of e via the relation (|1.6j) . This "reversion" can be quite nontrivial and is the subject 
of Section |1J Once this reversion is understood, one can then easily apply Theorem ll.2l 
to obtain the strong small-ball asymptotics. 

If Xfn is m-times integrated Brownian motion, then the lead-order asymptotics 
of — logP(y^ < s) have been studied by Chen and Li [J. In this one-dimensional 
case, lead-order asymptotics for the small-ball probability P(y^ < e) itself have been 
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studied by Gao et al. ([Sj, [H], and [7]), Nazarov [T^, and Nazarov and Nikitin jl4j. 
A closed-form expression for the Laplace transform of for any m has been derived 
in Gao et al. 0; we should also mention that in the special case m = 1 the Laplace 
transform of (|1.4|) had been obtained explicitly in Khoshnevisan and Shi ^U] several 
years earlier. 

We now turn our attention to higher dimensions. For d > 1, Csaki ^ computes 
the lead-order asymptotics of — logP(VQ^ < e), that is, the lead order for logarithmic 
small deviations of d-dimensional Brownian sheet. But a different approach is needed 
for nonzero m^'s. To illustrate, consider (mi, m2)-integrated Brownian sheet. Then 



oo oo 



h{x) = -logE exp <^ ^ ^ anun2^ni, 

[_ ni=l 712=1 



^ oo oo ^ 

(1-7) = Ei°Mi + 



2x 



ni=l ri2=l / z v / / 

where ani,n2 = 1/ [/^th ("^i)'^n2 ("^2)] and the Kn{m), n > 1, are the reciprocals of the 
eigenvalues of Am- The reason for moving to reciprocals is that the At's have nice 
representations as the zeros of certain entire functions If mi = 1712 = 0, then k„(0) = 
[(n — ^) vr] and Csaki's method works well because there are known expressions and 
bounds for the summations involved. However, if a component of m differs from then 
his approach is untenable. For example, if mi = and m2 = 1, then (k„(1)) are the 
solutions to the equation cos(z^/^) cosh(z-'^/^) + 1 = 0. When n is large, k„,(1) is quite 
close to [{n — ^) t^]^ (see, for instance, (51), but it does not seem possible to obtain 
explicit expressions. In this case we immediately see the complexity introduced when 
7712 = 1. 

Csaki handled only the case m = 0. Recently, Karol' et al. extended the classical 
result of Csaki ^ to obtain lead-order logarithmic small deviations for many Gaussian 
random fields; in particular, their treatment can handle our problem for arbitrary m. 
However, their methods do not seem to extend to obtaining asymptotic results beyond 
the lead order for the log small-ball probability. 

We note that the methods developed in this paper can be applied to many other 
Gaussian random fields besides m-integrated Brownian sheet X^. Consider a Gaussian 
random field whose covariance operator is the tensor product of marginal operators. 
If the Gaussian processes corresponding to these marginal operators each have a log 
Laplace transform whose asymptotic behavior can be established, then the methods 
of this paper will apply. Further, since (at least for X^) we can derive a complete 
asymptotic expansion for the log Laplace transform, we need only produce a complete 
asymptotic expansion extending Sytaja's theorem in order to obtain a complete asymp- 
totic expansion — not just the lead-order asymptotics — for the small-ball probability. We 
will produce the needed generalization of Sytaja's theorem, and apply it to small-ball 
probabilities for various Gaussian random fields, in future work. 

This paper is organized as follows. In Section[21we provide basic background concern- 
ing the Mellin transform. We derive asymptotic expansions for h{x) and its derivatives 
in Section 131 In Section 01 we discuss explicit expansions for the solution x* = x*(e) to 
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the equation h'{x) = e [recall (|1.6|) ]. Finally, in Section |21 we derive lead-order asymp- 
totics for the small-ball probability. 

Without loss of generality we will assume throughout this paper that 

< nil ^ ^2 < • • • < rrid < oo. 

For multibranched functions of a complex variable, we use the principal branch unless 
otherwise specified. 

2 Mellin transforms 

In this section we collect, for the reader's convenience, some useful facts concerning 
Mellin transforms. Our treatment follows closely the superb book manuscript by Flajolet 
and Sedgewick Hj. Throughout, we use the following notation for an open vertical strip 
in the complex plane: for real c, d, 

{c,d) := {s G C : c < Res < d}. 

2.1 Definition 

Suppose that / : [0, oo) — > C is locally integrable [i.e., is Lebesgue integrable over 
any bounded closed subinterval of (0, oo); this condition is met, for example, if / is 
continuous]. For those s E C such that x i— > f{x)x''~^ is Lebesgue integrable over 
(0,oo), we define the Mellin transform f* of / as the Lebesgue integral 

POO 

(2.1) r{s) := / f{x)x'^-Ux. 

Jo 

Keep in mind that Lebesgue integrability is a form of absolute integrability. 

2.2 Existence: the fundamental strip; transform of derivative 

Given local integrability, what is at issue for the existence of (|2.1() is the behavior of / 
near and near oo. It is easy to check that there exists a (possibly empty) maximal 
open vertical strip in which the integral 1)2.1(1 is well defined; that strip is called the 
fundamental strip. 

Suppose, for example, that / has the properties 

(2.2) f{x) = 0{x^) as X ^ 0+ and /(x) = 0(x") as x ^ oo 

with a < b. Let a denote the infimum of such a's, and P the supremum of such Vs. Then 
(— /3, —a) is a substrip of the fundamental strip [and equals the fundamental strip in 
typical cases, such as when f{x) = (l + o(l))cox^ as x ^ 0+ and f{x) = {l + o{l))cooX°' 
as X — > oo with co,Coo ^ 0]. If, further, / is continuously differentiable and monotone, 
then one can check using integration by parts that {1 — (5,1 — a) is a substrip of the 
fundamental strip for /' and that 



(2.3) 



if'Tis) = -is - l)f*{s - 1), s€{l-P,l-a). 
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In this paper, we will make key use of the following Mellin transform pair: 

vr 



(2.4) /(a;) =iog(l + x) ^r(s) 



s sin(7rs) 



with fundamental strip (—1,0). For justification of the double arrow, see the start of 
the next subsection. 

2.3 Inversion and the mapping property 

Suppose for simplicity that / is continuous and has a nonempty fundamental strip, 
which we denote by (— /?, —a). Then for any c in the interval (— /3, —a) we have 

1 pc+ioo 

(2.5) /(:c) = — / ns)x-'ds. 

This establishes a correspondence between functions and their Mellin transforms. 

Mellin transforms are useful in asymptotic analysis because, under a certain growth 
condition on the meromorphic continuation of the Mellin transform /*, there is a cor- 
respondence between the asymptotic expansion of the function / near or oo and the 
singularities of f*{s) for s G C. (For our purposes, we need only asymptotic expan- 
sions near co, so those are all we will discuss here.) Flajolet and Sedgewick [4^ call this 
correspondence the mapping property of Mellin transforms. To explain more fully, it is 
convenient to introduce the following formal-sum shorthand notation. 

Let (/? : — > C be a meromorphic function in a domain 17 and let C 17 be the set 
of its poles. For each sq G V, let Asq(s) be a truncation of the Laurent series for ip{s) 
at So S V] we assume that As(,(s) contains at least all terms of order (s — sq)~^ with 
A; > 0. We then write 

(2.6) <^(^)^ E [^«o(^)]s=so (se^^), 

or, when no confusion results, simply 

^{s) X ^so{s) {s G fi), 

and call (|2.6j) a singular expansion of 93 in J7. 

The following two theorems are essential for our results. The first of these theo- 
rems gives (among other things) sufficient conditions for meromorphic continuation of 
the Mellin transform /*. The second theorem shows that a singular expansion of /* 
characterizes the asymptotic behavior of /. 

Theorem 2.1 (Direct Mapping Theorem). Let f have Mellin transform f* at least 
in a nonempty strip {—13, —a). Assume that f{x) admits as x —> 00 a finite asymptotic 
expansion of the form 

(2.7) f{x)= c^,kxH^ogx)'' + 0{x^), 
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for a finite set A of pairs {(,,k), where the ^'s satisfy 7 < ^ < a and the k's are 
nonnegative integers. Then f* is continuable to a meromorphic function in the strip 
{—(3, —7), where it admits the singular expansion 

Theorem 2.2 (Reverse Mapping Theorem). Let f be continuous in (0, cxd) with 
Mellin transform f* existing at least in a nonempty strip (— /?, —a). Assume that f*{s) 
admits a meromorphic continuation to the strip (— /?, —7) for some 7 < a with a finite 
number of poles there, and is analytic on Re s = —7. Assume also that there exists a 
real number rj E (a, 13) such that 

(2.8) f*{s) =0{\s\-'^), with r>l, 

when \s\ — > 00 in —rj < Res < —7. ///* admits the singular expansion 

with all k's in the sum nonnegative integers, then an asymptotic expansion of f{x) as 
X ^ 00 is 

f{x)= E c^,kX^{logxf + 0{x'). □ 

2.4 Termwise differentiation of an asymptotic expansion 

The proof of the foUowing resuh is a simple exercise in the use of the two mapping 
theorems. Note the condition r > 2. 

Corollary 2.3. Let f be continuously differentiable and monotone, and let a and (3 
be defined as in the paragraph containing 1)2. 2() and (|2.3() . Assume that a < (3, that 
f{x) admits as x ^ 00 a finite asymptotic expansion ()2.7|) as in the statement of the 
Direct Mapping Theorem, and that f* has a meromorphic continuation that is analytic 
on the line Res = —7 and satisfies the growth condition 1)2. 8|) with r > 2. Then the 
expansion (|2.7|) can be differentiated termwise: as x —> 00, 

f{x)= E c^,kX^-\logx)^-\Clogx + k) + 0{x^'^). □ 

2.5 Harmonic sums and the separation property 

The mapping property discussed in Section 12., SI is particularly effective for asymptotic 
analysis of a harmonic sum 

(2.9) F(x)=EW(^fcx) 

k 
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with /ifc > for every k. The reason is the readily checked separation property that the 
Mehin transform of (|2.9|) has the simple product form 

(2.10) F*{s) = ns)Y,Xkf^^\ 

k 

valid for s in the intersection of the fundamental strip of / and the domain of absolute 
convergence of the generalized Dirichlet series A^^^**. 



3 Asymptotic expansion for the log Laplace transform 

Recall the expression H1.7() for the log Laplace transform = h oi V'^. To avoid 
nuisance factors it is convenient for us to reparameterize hm by defining, for x > 0, 

(3.1) -21„gE exp {-|^^} ^ " E (' - ^ 

Notice that h^{x) = ^Lin(23;). In this section we derive a complete asymptotic expan- 
sion for Lm, for general m. First we treat the one-dimensional case. Then we bring 
Mellin asymptotic summation to bear. That is, we factor the Mellin transform in 
terms of the transforms L*^^^ using the separation property, analyze the singularities of 
the factors using (among other things) the Direct Mapping Theorem, and finally obtain 
the asymptotic expansion for Lm (or for h^) using the Reverse Mapping Theorem. Us- 
ing Corollarv 12.31 we also derive corresponding expansions for the first two derivatives 
of /im, as needed for our application in Section [31 of Sytaja's Tauberian Theorem. 



3.1 The one-dimensional case 

We begin by analyzing m-integrated Brownian motion, i.e., the d = 1 case 



(3.2) 



CO 



n=l 



log 1 + ^- 



For integer i, < i < 2m + 1, set ivi := exp{i7r^^^} = loi, vg := exp{i7r 2m+2 i ~ 
VQUjg, and I3g{x) = /?£ := x^/'^^^+^^zv^. We employ the following result from Gao et al. 
(cf. Theorem 6 in f^), which gives an exact expression for Lm{x): 

Lemma 3.1. For x > we have Lm{x) = log | det N{x)\ — (m + 1) log(2m + 2), where 



( 1 

(Jo 



N(x) :-- 



1 



\ ^0 



1 \ 



m+1 /32m + l 

'.^2m+l^ 



2m+l /32,„+i / 
'^2m+l^ / 



□ 
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The following important asymptotic consequence will be needed. 
Lemma 3.2. As x ^ oo, the quantity Lm{x) defined at (|3.2|) has the expansion 



CSC 



vr 



„l/{2m+2) 



2m + 2 
+ O ( exp <j — sin 



+ 2 log I det U\-{m + 1) log(2m + 2) 



vr 



2m + 2 

where Um = U is the Vandermonde matrix 



X 



l/{2m+2) 



(3.3) 



U :-- 



/ 1 1 



\ ^0 '^l 



1 

















, ,m , 



Proof. We need to study the large-x behavior of N{x). Multiply the last m + 1 columns 
of N{x) by e^°, e^^, e^™, respectively, and use Pj = —Pm+i+j to obtain the following 
matrix: 



/ 1 



N{x) 



^0 

LO^+^ef^ 

For < £ < m we have 
and therefore 



1 

m+l f3„ 



2m+l (3„ 



m+l 



m+l 



1^, 



2m+l 
m+1 



g/3m 

i^2m+ie' 



\ 



^2m+l^ 
m+1 



UJ. 



2m+l 



, ,2m+l 
"^2m+l 



exp{— sin {-n /{2m + 2)) x 



l/(2m+2) 



} 



detiV(a;) = det iV(oo) + 0(|e^° |), 



where N{oo) is the matrix obtained from N{x) by replacing each entry with a factor 
e^* by 0: 



iV(oo) 



/ 1 • 


• 1 





• 









• 





■ 




m+1 
■^m+l 


• 

m+1 
"^2m+l 


Vo • 


• 


, ,2m+l 
'^m+l 


, ,2m+l 
■ "^2m+l 



A simple calculation reveals that | det A^(cxd)| = | det C/p 7^ 0, where U is defined at 1)3. 3(1 . 
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Finally, using 



g/30 g/3l . . . g/3rj 



we obtain 



Lmix) + (m + 1) log(2m + 2) 



exp < — CSC 



vr 



2m + 2 



X 



l/(2m+2) 



log I e 



-/3o^-/3i 



• • e 



log 



CSC 



CSC 



exp CSC 
vr 

2m + 2 
vr 

2m + 2 



vr 



^™ det A^(x)| 

:e1/(2™+2) UetiV(x) 



2m + 2 

xy(^^+^) + log\detN{oo)\+0{\e^°\) 
^i/{2m+2) + 2 log I det U\ + 0( exp | - sin 



vr 



2m + 2 



X 



l/{2m+2) 



as desired. 



□ 



Remark 3.3. From Lemma 13.21 we see that L^(x)/x^/(2m+2) ^ f^j^-^g 

nonzero limit 

as X — > oo. Also, it is easy to see that Lm{x)/x has a finite nonzero limit x 0+. 
Therefore, the fundamental strip for the Mellin transform of Lm is (—1, — l/(2m+2)). 

Remark 3.4. As a consequence of the Direct Mapping Theorem we see that is 
continuable to a meromorphic function in the strip (— l,oo), with 



(3.4) Ll,{s) 



CSC 



2m+2 



S + 



1 



2m+2 



+ 



-2 log I det U\ + (m + 1) log(2m + 2) 



s=0 



in particular, in this strip L^{s) has only simple poles, at s = — l/(2m + 2) and at 
s = 0. 

Remark 3.5. Although our treatment keeps m fixed, the reader may be curious as 
to the large-m behavior of the term 2 log | det Um\ appearing in Lemma 13.21 Using the 
well-known formula for the determinant of a Vandermonde matrix and approximating 
the double sum that results after taking logarithms by a double integral, we find 



m 



1 



■ log 2 + 



2 log I det Ur, 
We omit the details. 

c 

Since Lm{x) = ^ log ( 1 + 



0<x<y<l 



vr 



log sin \ —{y — x)] dxdy 



7C(3 ) 
2vr2 



-m 



n=l 



(3.5) 



K„(m, 

-^m(-s) = 



, it follows from (fTl)) and (fTTU)) that 



vr 



ssin(vrs) ^ 

n=l 
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or what is the same, 

oo 

(3.6) KM-=^<(.rn) 



n=l 



S Sm TTS 



vr 



Since, as discussed in the next subsection, Kn{m) grows hke n^™"^^, (|3.5j) - (|3.6j> are vahd 
(without continuation) in the strip (—1, — l/(2m + 2)). 

3.2 The generalized Dirichlet series Km{s) 

When we apply the Reverse Mapping Theorem to in Section 13.31 we wih need to 
verify the growth condition ()2.8|) . This will rely on corresponding estimates for the 
one-dimensional LJj^, and so we need to study the growth of L'^{s), or equivalently of 
the generalized Dirichlet series Km{s). 

We can do this using the following very sharp asymptotic estimate for Kn{m) (see 
also Gao et al. |Sj). 

Lemma 3.6 (Gao et al. Theorem 2). For m = we have Kn(0) = [('^ ~ ^) 
for every n. For each fixed m > 1, as n ^ oo we have 



K„(m) = [(n - i) vr] 



2m+2 



1 + 



n ^ exp 



-vr sm 



vr 



m + 1 



n 



□ 



From Lemma 13.61 it is easily seen that the domain of absolute convergence for Kn 
is the strip {—oo, — l/(2m + 2)). In this strip we can write 



Kim)] 



where 



m) 



n=l 



n=l 



( 7r \ (2m+2)s 
k2) 



vr 



(2m-|-2)s 



C(-(2m + 2)s) 



with kn{m) '■= [{it- ~ ^) T^] i here (" denotes Riemann's zeta function. Using Lemma 
13.61 it is easy to check that the remainder series Yl'^=i ['^n("^) ~ '^n("^)] converges abso- 
lutely for all s € C and defines an entire function of s that is 0(|s|) in any strip {—R, R) 
with < -R < oo. The generalized Dirichlet series Km{s) is meromorphically continu- 
able for s S C, with a single simple pole at s = — l/(2m + 2). In any strip {—R, R), the 
continued Km{s) grows at most polynomially in |s| (see, e.g., Section 5.1). Putting 
the pieces of the above argument together, we have established the following result. 

Lemma 3.7. For any < R < oo, the meromorphic continuation of the generalized 
Dirichlet series Km{s) at (|3.6() grows at most polynomially in \s\ as \s\ — > oo in the 
strip {-R,R). □ 
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3.3 The general d-dimensional case: separation 

With results regarding the one-dimensional processes in hand we now turn our atten- 
tion to arbitrary dimension d and vector m. The following simple lemma exploits the 
harmonic-sum structure (|3.1|) of I/^(s). 

Lemma 3.8. For any dimension d>l, 

(3.7) ^.^(,)^(ii^Y"{lLl,,(s), 

valid (without continuation) in the strip (— 1,— l/(2mi +2)). 

Proof. When d = \ the result is trivial (recalling Remark 13. 3|) . For d > 2, we observe 



Lm{x) — V] ^ ( I 



s sin(7rs) 



vr 



whence, by the separation property (|2.1U|) . (|3.6() . and induction, 

oo ' ( \ 

nd=l 

d~l d 

i=i 

for s € (—1, — l/(2mi + 2)), as claimed. □ 

In preparation for our main theorem (Theorem I3.1U|) we next consider the singular 
expansion of the meromorphically continued L'^{s). In light of Lemma 13.81 it is enough 
to understand the meromorphically continued L^.{s) for each j. Each L^.(s) has a 
simple pole at s = 0, but ssin(7rs) has a double zero at s = 0; thus, if d > 2, then L^{s) 
has as singularities in (—1, oo) only poles at s = — l/(2mj + 2) for j = 1, . . . , d. 

Let 



CSC ( — - — ] x^/(2m+2) + 2 log I det U\ - (m + 1) log(2m + 2) 
2m + 2 



denote the exponentially small remainder term in the asymptotic expansion for Lmix) 
in Lemma 13.21 From Lemma 13.21 and (the proof of) the Direct Mapping Theorem, for 
s / — l/(2mj + 2) we have [compare (|3.4j) ] 



CSC 



(3.8) r^^is) = ^'T'' +Am,{s). 

Here 

-2 log I det Umj I + {ruj + 1) log(2mj + 2) 



(3.9) A^^ (s) = [ Lm^ {x)x'-^ dx + 
Jo 



f oo 

+ / rm^{x)x''~^ dx 



is analytic in the strip (— 1, cxd) except at s = 0. 
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Remark 3.9. From the above discussion and Lemma 13.81 we immediately conclude 
for d > 2 that L'^{s) is meromorphic in ( — 1,cxd) with its only poles located at s = 
— l/(2mj + 2) for j = 1, . . . ,d. Further, Lemmas 13.71 and 13.81 together with the expo- 
nential growth of the sine function along vertical lines in C, imply that L'^(s) satisfies 
the growth condition (|2.8|) for every 1 < r < oo. 

3.4 Main theorem: Complete asymptotic expansion of L^iix) 

With the previous results we are now able to obtain a complete asymptotic expansion 
for Lm{x). In order to state our main result, we suppose 

mi = ■ ■ ■ = < mt-^+i = ■■■ mt^+tz < ■ ■ < mt^+...+tg-i+i = ■ ■■ = mtj+...+tg, 

where the number of groups of ties is g > 1 and, for 1 < i' < g, there is a tie of size 
> 1 in group v; thus d = ti + ■ ■ ■ + tg. We will denote the common m-value in the 
vth group by fh^. Set 

. - 1 
■ 2m^ + 2 ' 

In order to use the Reverse Mapping Theorem, we need the singular expansion of L^{s) 
near each pole s = {u = 1, . . . , g). From (|3.7|) and (|3.8|) . near s = —^y we have 

(3.10) lAjA = (ifi^M)"' j|n [ijj.)]'") 

The product of the first two of the three factors on the right in ()3.1U|) is analytic for s near 
and therefore has a Taylor expansion that can be computed up through the term 
of order (s + ^jy)*""^; in regard to this computation, note that repeated differentiation 
under the integrals in (|3.9|) is easily justified. Similarly, we can expand the analytic 
function Ain^{s) up through order (s + ^uY"~'^ and use a mulitnomial expansion to 
obtain a Laurent series for the third factor up through order (s + ■ Multiplying 
these expansions together we can get a Laurent series for (|3.1Up up through the term 
of order (s + ■ Applying the Reverse Mapping Theorem (the growth condition is 
satisfied: see Remark 13. 9|) we arrive at our main result. 

Theorem 3.10. For any m-integrated Brownian sheet in dimension d > 2, Lm(3^) has 
as X ^ oo a complete asymptotic expansion of the following form: 

g U-i 

(3.11) Lm(x) = ^Y1 c.,fcX«^(logx)^' + 0(x-^) 

for any R > 0. The values ofc^^k ore computed as outlined in the previous paragraph. □ 

Remark 3.11. In fact, the error term in Theorem I3.1UI is exponentially small in a 
positive power of x. We have written a complete proof, but will provide only a sketch in 
the next paragraph (as the details are straightforward but laborious), that the remainder 
term is bounded by exp i^—cx^/^'^^'^'^d+'i)) | fQj- g, certain constant c depending on (d and) 
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m. We have not tried to optimize this bound; in particular, it may be that the power 
l/{d{2md + 2)) can be improved to l/[2(mi + • • • + m^^) + 2d]. 

The Reverse Mapping Theorem 12.21 is proved by invoking the inverse Mehin trans- 
form formula ()2.5() with c = —r] and then shifting the line of integration rightward to 
Re s = —7 by means of the residue theorem; the growth condition (|2.8j) is used to justify 
this shift rigorously. One then uses the growth condition again to bound the shifted 
integral by 0{x'^). When sharper growth estimates of f*{s) are available, as they are 
in our case / = Lm (cf. Remark 13. 9() . one can let I7I grow with x and obtain sharper 
remainder bounds. Our proof for Lm uses well-known growth estimates for the Riemann 
zeta function in vertical strips and takes 7(0;) = — cx^/^'^*^^'"'^"'"^-*) -|- 0(1) for a suitable 
positive constant c. 

Because h{x) = ^Lm(23;), it follows from Theorem I3.1UI that h[x) also has an ex- 
pansion of the form ()3.11|) . namely, for any R > 0, 

(3.12) h{x) =^211 ^-.'=(0) ^^"(logx)^ + 0(x-^) 

u=l k=0 

where 

c,,fc(0) :=2«''-i J^c,,,Q(log2)^-*^ 

(and error term as refined in Remark I3.11|) . In order to apply Sytaja's theorem in 
Sectional we will need the following extension of Theorem 13. 101 to the derivatives of h. 

Lemma 3.12. For j = 0, 1, 2, ... , the function h^^^ has as x —>■ 00 a complete asymptotic 
expansion of the form 

9 U-i 

(3.13) h^^\x) = ^Y1 c^,fc(j>^''"^'(logx)'= + 0(x-(^+^)) 

j/=l k=0 

for any < i? < 00. The expansions are obtained by successive termwise differentia- 
tions, starting with ()3.12l) . In particular, with Cy^kij) '■= whenever k > ty, 

Cz/,fc(l) = ivCi,,k{^) + {k + l)c,,,fc+i(0), 
c^,fc(2) = {iu - l)c^,fc(l) + {k + l)c^,fc+i(l) 

= Uiu - l)c.,fc(0) + (2^, -l){k + l)c,,fc+i(0) + {k + 2){k + l)c,,fc+2(0), 

and, for any < R < 00, 

a t^-i 

h{x) - xh'ix) = J] J] [(1 - 6)c.,fc(0) -{k + l)c,,fe+i(0)]x«''(logx)^ + 0(x-«). 

u=l k=0 

Proof. Denote the eigenvalues of Am by a^, as at (|1.5() . Then h{x) = log(l-l-anx) for 
all X > 0, and the dominated convergence theorem justifies the termwise differentiation 
giving, for j = 1, 2, . . . and x > 0, 

h^^\x) = i-iy^Hj - 1)! E «n(l + «nX)-^\ 
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In particular, each function h^^^ is monotone; and as x — > 0+ we have the simple 
estimates 

h{x) = (1 + o(l)) X a„, h^^\x) = (1 + o(l))(-l)^-i(j - 1)! <■ 

n n 

It is then easy to use induction on j, together with Corollary 12.31 and Remark 13.91 to 
complete the proof. □ 

Remark 3.13. For each j = 0, 1, 2, . . . , the error term in (|3.13|) is exponentially small 
in a positive power of x; compare Remark 13.111 

3.5 Expansion of Lm{x): examples 

We next give three simple examples of the computations entering into the expansion in 
Theorem 

3.5.1 Lead-order asymptotics in the general d-dimensional case (d > 2) 

Suppose d>2 and 

m := nil = ■ ■ ■ = rrit < mt+i < • • • < ma- 
According to Theorem 13. lUl the lead-order asymptotics for Lm{x) are given by 

Lm{x) = (l+o(l))ci,t.ixl/(2™+2)(logx)*-l, 

since (as will become evident below) ci^t~i > 0. This lead-order term corresponds via 
the Reverse Mapping Theorem to the term of order (s-|-^i)~* in the expansion of L'^{s) 
near its pole at s = — ^i. From (|3.1U|) we see immediately that, as s — > — ^i, 

'6sin(7r^i)\'^"W -rr r* / ^nI /-csc«i)^* 



^;.(^)=(i+o(i))(^i^^)" n ^^.(-^1 



s + i 



j=t+l 

d-l-t 

1 \ \ (-l)*-i(t-l)! 



^ + [(2m + 2)vr]^-i(t - 1)! 11^ "^^^ 2m + 2)] {s + 6) 



j=t+i 



from which it follows that 



(^•1^) -^---^= [(2m + 2)vr/l(t-l)! n ^ 



1 



[(2m-F2)7r]^-i(t- 1)! 1 .jLi ™j 2m + 2 

Thus, in general, computation of the lead-order term for L^{x) requires (numerical) 
evaluation of the fundamental-strip Mellin transform values 

(3.15) Ll,^ [-^^^ = ["^-^ ^-(2n.+3)/(2„.+2) > o_ 
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However, in the particular equal-m's case where t = d we obtain the simpler result 

L (x) - (l + o(l)) ^^^(2;^) x^/(2m+2)Q xd-l 

Remark 3.14. From (|3.6jl we see that an alternative to numerical evaluation of the 
integral (|3.15|) is numerical evaluation of the generalized Dirichlet series (^~ 2m+2 ^ ■ 
Evaluation of such sums arises in the approach of [2], which is more computationally 
intensive since it involves numerical root-finding. 

3.5.2 Complete asymptotic expansion: distinct m^'s 

Given Theorem l3.1Ul it is almost trivial to obtain a full asymptotic expansion for Lmix) 
when the vector m has distinct components. That is, suppose that g = d > 2, so that 
tiy = 1 for 1/ = 1, . . . ,d. Then Theorem 13. lUI implies 

d 

(3.16) L^{x) = Y,Cu xi/(2™.+2) ^ o{x-^) 
for any < i? < cxd, where 

r / \ld-2 

(3.17) C, 



[{2m^ + 2)7t' 



\d-l 



3.5.3 Full asymptotic expansion: d = 2 

Theorem 13.101 provides a complete asymptotic expansion of LIy^{x) for any vector m. 
However, calculations quickly become quite cumbersome to perform by hand; symbolic- 
manipulation software such as Mathematica or Maple can then be of great help. Com- 
plementing the above result for distinct m's, here is another case where we can spell 
out the expansion without too much notation. 

Suppose that m := mi = m2, i.e., that g = 1 and d = ti = 2. We then find 

Lm{x) = Cl,l xl/(2™+2) log 3. + ci,o X^l^^'^+^'^ + ©(x"^) 

for any < i? < cxd, where 
(3.18) ci,i 



(2m + 2)7r 



2m+2) ^^^(2m+2j ( 2-m+2 j \^2m+2 
(3.19) Cifi = — \ \ — , 

[m + Ijvr TT 2m + 2 

and Am{-) is computed via (|3.9() . When m = in this example, i.e., when V is the 
L^-norm of two-dimensional Brownian sheet, we have from Lemma l3. II that 



Lq{x) = logcosh(a; 



l/2x 
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and 



Ao{-l/2) = j Lo(x)x"^/2dx-21og2+y Lo(x) - x^/^ + log 2 x~^''^ dx 



-0.3624. 



Thus, 



^(0,0) {x) = log X + ci,o x^/" + 0( 



1/2 , 



27r 



for any < i? < oo, and here 



l + ^o(-l/2) 
ci,o = — ^ 0.2029. 



vr 



4 Reversion of h'{x) = £ 

In principle we are now in position to use Sytaja's theorem in conjunction with the 
asymptotic expansions for h{x) — xh'{x) and h"{x) in Lemma [3.121 to state a strong 
small-deviations result for integrated Brownian sheets. The problem is that, for more 
explicit results, we still must solve h'(x) = e to get x* = x*(e) in order to obtain the 
needed expansions in e, up to an additive term o(l) for h{x*) — ex* and up to a factor 
l + o(l) for h"{x*). 

In this Section 0] we discuss the reversion of the asymptotic expansion for h'[x) 
given in Lemma 13.121 We begin in Section f4. II bv showing that the exponentially small 
remainders in all our asymptotic expansions may safely be ignored when carrying out 
reversion and applying Sytaja's theorem. In Section 14.21 we determine the lead-order 
asymptotics for x*. Clearly, the lead-order asymptotics agree for x* and the (unique 
real) solution xq = xo(e) to the equation 

ti-i 

(4.1) e= 5^01,^(1) x«i-i(logx)'= 

fe=o 

obtained from the equation e = h'[x) by including only those terms with v = 1 (corre- 
sponding to the largest power of x) in the expression 1)3. 13(1 (with j = 1) for h'[x). Next, 
in Section U31 we obtain a complete asymptotic expansion for x* in terms of elementary 
functions and xq. Finally, in Section f4. 41 we discuss exact computation of xq. 



4.1 Truncating the asymptotic expansion for h' suffices 

Our next result implies that, as far as obtaining lead-order asymptotics for the small- 
ball probability, we may act as if the asymptotic expansions in Lemma 13.121 were exact 
expressions. 

Lemma 4.1. Let h denote the truncated asymptotic expansion for h: 

g tu-i 

i/=l A:=0 
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and let x satisfy h[x) = e. As e ^ 0+, the following errors all tend to zero faster than 
any power of e: 

x — x*, h^^\x) — h'^^\x*) (for j = 0,1, ... ), [h{x) — xh' (x)] — [h{x*) — e x*]. 

Observe that h'{x) decreases monotonically to for sufficiently large x; thus for 
sufficiently small e there exists a unique solution x = x{e) to h'{x) = e. We remark in 
passing that h^^^ can also be obtained by truncating the asymptotic expansion for h^^'^ 
given in Lemma 13.121 

Proof. From definitions and Lemma l3.12l we have 

h'{x) = e = h'{x*) = h'{x*) + 0{{x*)-^) for any < < oo, 

which by the simple Lemma 14.31 below can be written 

h'{x) - h'{x*) = 0{e^) for any < i? < cx), 

as e — > 0+. It follows using the mean value theorem that 

X — X* = 0{e^) for any < i? < oo. 

The other assertions follow readily. □ 

Remark 4.2. The errors in Lemma [4.1l are each exponentially small in a positive power 
of 1/e. This assertion follows easily from the case j = 1 of Remark 13.131 

4.2 Lead-order reversion in the general rf-dimensional case (d > 2) 

If all we want is a standard "weak" small-deviations result, we need only obtain the 
lead-order asymptotics for x*{e) [equivalently, for xo{e)]. This is easy. Recall from 
Lemma r3. 121 that h'{x) = (1 -|-o(l)) ci^f_i(l) x^~^ (log where t := ti is the number 

of mj's equal to the smallest value m := mi and ■= (,i = l/(2m -|- 2). The coeffi- 
cient ci^t_i(l) is given by ci^t_i(l) = ^2^^^ci^t-i, with ci^t-i given by 1)3. 14|) . The proof 
of the following lemma is left to the reader. 

Lemma 4.3. // h'{x*) = e and xq is the solution to (|4.1j) . then as e ^ 0+ we have 
X* = {1 + o(l)) xo = (1 + o(l)) Xq, where 



(4.2) xo = ioie) 



ci.t-i(l) 1 A 1^*"'' 



1 



□ 



We will use Lemma f4.31 in Section|Slto obtain a rather simple "weak" small-deviations 
result for arbitrary m; see Theorem 15.11 Lemma 14.31 also provides us with all the 
asymptotic information about h"{x*) we need in applying Sytaja's theorem: 

Corollary 4.4. If h'{x*) = e, then, with xq defined at (|4.2j) . as e — > 0-|- we have 



-{x*fh"{x*) = {l + o{l)){l-i)x,e 
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Proof. This is routine. Recall Lemma l3 . 1 2 1 and in particular the coefficients 
ci,t_i(2) = e(C - l)2«-ici,t_i = -(1 - Oci,t-i(l). 

Then 

-{x*fh"{x*) = -{l + o{l))4h"{xo) = -(l + o(l))ci,t_i(2)i« (logfo)*-' 
= -(l + o(l))ci,t_i(2)4 (^^logi) , 

and it is easy to check from definitions that 

-ci,t_i(2)4 (^^log^^ ={l-^)xoe. □ 

4.3 Reversion: an asymptotic expansion for x* 

We next derive a complete asymptotic expansion for x* in terms of elementary functions 
and the solution xq to H4.1|) . We will discuss exact computation of xq in Section 

Dropping the error term from the expansion (||-{.ll-i|) (with j = 1) for h'{x) (as justified 
by Lemma l4.1() . let us write the equation e = h'{x) in the form 

(4.3) e = 5]/,(x), 
where 

(4.4) /^(x) :=x-'''' ^ a^,fc(loga;)^ u=l,...,g; 

here, for abbreviation, we have set ^ := ^ (1) andr/^ := 1-^^ = (2m^ + l)/(2m^+2), 
so that 

1/2 < r?i < r?2 <•••%< 1. 

(The notation rhi, is as in Section 13.41 ) g = 1, then we have simply x* = xq; so we 
assume now that g > 2. 

The main result of this subsection is the following complete asymptotic expansion 
for X* /xq in terms of inductively defined quantities yj (j = 0, 1, 2, ... ). Let yo := 1- 
Suppose that j > 1 and that we have defined yo, . . . , yj-i- Then set 

Vj'-i ■=yo + yi-\ h Vj-i and := xq y^_^ 

and define 

yt-i [Et=iMxj-i)-e] 



1) 
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Proposition 4.5. For each j = 0, 1, . . . we have 

(4.6) — = 1 + yi + ya + • • • + (1 + o{l))yj 

xo 

and 



(4.7) 



yj = 0\ [en ^ ( log 




where A := max{U - I - {h - 1)^ ■■2<u<g}. 

Note that, with x_i := 0, (|4.(ij) can be written equivalently as 

(4.8) X* = x,_i + (1 + o(l))xoy,- (j = 0,l,...). 

Before we prove the proposition, we illustrate its apphcation in two special cases. 

Example 4.6 (ti = 1). Suppose ti = 1, whence /i(x) reduces to aifiX^"^'^ and it is 
elementary that xq = (aifi/e)^/^''- . In this case, (|4.5|) simplifies to 



(4.9) 



y 



m«i,o 



1 -e 



u=l 



If we assume further, as in Section [3. 5. 21 that all the m^'s are distinct (i.e., g = d), then 
fuix) simplifies to a^fiX'^" and (|4.9|1 can be written in the form 



(4.10) 



^a.,oa-2^/^^e(''^/''^)-i(y+_i)-''^-l 



.u=l 



From this it is easy to prove by induction that each yj has an asymptotic expansion in 
increasing powers of e wherein each power is a nonnegative integer combination of the 
numbers (rju/Tji) ~ 1 with 2 < v < d. The same is therefore true of x* /xq. 

Example 4.7 (mi = 0). Suppose mi = fhi = 0, whence ryi = 1/2; since fh2 > 1, we 
also have r/2 > 3/4. In this case it is not hard to see that the finite expansion 



- = 1 + yi + y2 + 0(^3) = 1 + 2/1 + y2 + of e^''^-^ fbg - 

Xo \ \ e 



i + yi + y2 + 0^e'/' (log^) ^ 



is sufficient to obtain the lead-order asymptotics for the small-ball probability. To be 
concrete, suppose further that there are no ties, so that we are in the context of 1)4. 10() 



20 



and xq = (aifi/e)'^. Then 

d 



u=2 



yi = 2^a.,oaJ^^e2,.-i^ 



y2 



y=2 



u=2 



2r?„ J2r)y-1 



6r;2-3\ 



(4.11) 



where the first remainder term in (|4.11|) is used if d > 3, and the second if d = 2. So, 
whether ri2 > 3/4 or 7/2 = 3/4, and whether d > 3 or d = 2, 



Xo 



l+yi + o{e) = 1 + 2 ^ a,,o a J'^'^e^r,.-^ + o{s), 



u=2 



which we wiU see in Example 15.81 is sufficient to obtain the lead-order asymptotics for 
the small-ball probability. 



Proof of Proposition \4.5{ We prove (|4.8jl and (|4.7)) together by induction on j. The 
base case j = of the induction is simple. For j > 1, our induction hypothesis is that 
the {j — l)st instances of (|4.8j) and (|4.7|) hold; in particular, we can write 



X* = Xj_2 + xoyj^i + xoy = Xj-i + xoy = xj-i 1 + 



y.i-i 



for some y = o(?/j_i). For = 1, . . . , g' we then have, from 1)4.4(1 . 



f4x*) = fJxj-i 1 + 



yj-i. 



1 + 



k=0 



logx,_i -I- log 1 -I- 



For > 2 we conclude 



f,{x*) - = 0(y/,(x,_i)) = 0(y V^(logxo)*''-^) 

(4.12) = O U e''''/'?M log - j \=0\y e'^^/^i / log - 

where the third equality follows from Lemma 14.31 For = 1 we find 
(4.13) 

y 



/i(^*)-/i(a;,-i) = -(i + o(i))- 



y 



ti-i 



m/i(2;i-i) - ^ A;ai,fc(logXj„i)'' ^ 



fc=o 
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which we note is of the same order as ye. Therefore, 

9 9 

e = ^Mx*) = ^/,(x,_i)-(l+o(l))^ 
u=i u=i 



ti-i 



"j-i) 



fc=0 



and so [recaUing (|4.5p ] y = (1 + o{l))yj. This estabhshes the jth instance of ()4.8() . 

To establish the jth instance of H4.7() , we note again that the denominator of (|4.5|) is 
precisely of order e. We must also estimate the difference appearing in the numerator. 
Suppose that j > 2; the estimation for j = 1 is similar and left to the reader. By slight 
modifications of the arguments for ()4.12() - H4.13() . for u >2 we have 



(4.14) 

and, also utilizing ()4.5() . for u = \ we have 



/l(Xj_l) - /l(Xj_2) = 

9 



ti-1 



fl{Xj-2) - XjH^ X] ^«l,fc(logXj_2)'' ^ 



fe=0 



(4.15) 



Xj-2) - e 



+ 0(2/1-1^) 



Summing the g equations (|4.14|) - H4.15() and rearranging. 

By the induction hypothesis and our assumption that j > 2, the first of the two O(-) 
terms predominates and the estimate 



is established, completing the induction. 
4.4 Computation of xq 

There remains the task of solving the following equation to obtain xq: 



□ 



(4.16) 



e = h{x) = ^flfex ''(logx)*^, 



fc=0 



with t := ti, ak := ai^k = ci,fc(l), and r? := r/i = 1 - ^1 = (2mi + l)/(2mi + 2). As 
has already been noted in Example 14.61 this is trivial if t = 1: then xq = {ao/e)^^'^. 
For general t, of course, one can resort (for given e > 0) to Newton's method or other 
root-finding methods. What we will derive in this subsection is a series representation 
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of the solution, from which wih fohow an asymptotic expansion for xq sufficient to yield 
a complete asymptotic expansion for the logarithm of the small-ball probability, but not 
for that probability itself (for which the exact value of xq should be used). 

It is instructive to begin by considering the case t = 2. Setting w := —r][logx + 
(ao/ai)], equation (|4.16|) can be rewritten as 

we^ = exp{— r/ao/ai}e. 

ai 

This equation for w cannot be solved in terms of elementary functions, but it does have 
the solution 



w 



W 



Tj 

— exp{-r/ao/ai}e 
ai 



in terms of the Lambert T^-function. (See (2|; more specifically, W here is the branch 
W-\ in the notation there. We could alternatively work in terms of the closely related 
glog function defined in jH].) According to the paragraph following equation (4.20) in [2] 
we have, in terms of z := (ry/ai) exp{— r/ao/ai}e, the expansion 



(4.17) 



w = Wi-z) 



^ CO OO y 1 \ ^ 

log Z - log log - + ^ ^ drs I log log - j (log . 



\-(r+s) 



in [IH) by drs = hi-iy 



As remarked in j2|, the series in (|4.17|) is absolutely 



with drs given in terms of Stirling numbers of the first kind (see equation 1.2.9-(26) 

r + s 
r + 1 

and uniformly convergent for small e and also serves as an asymptotic expansion for 
w — log z + log log i as e ^ 0+. It is not hard to show that this result can be rearranged 
to one of the following form, wherein the series has the same convergence and asymptotic 
expansion properties as in 1)4. 17(1 : 



(4.18) xo = xo X 
Recall that, in general, 

(4.19) 



OO CO 



l + ^2^dr 

r=0 s=0 



log log 



at-1 1 



V 



t-i 



- • log - 



log 



i-1 



-(r+s)' 



with di 



00 



0. 



1/77 



is the lead-order approximant to xq found in Lemma 14.31 and for t = 2 this reduces 



to Xq 



21 . 1 . log 1 



i/'n 



. In particular, the first correction term in (|4.18j) then has 

coefficient doi = doi/r] = l/rj. 

The expansion 1)4. 18(1 can be extended to general values of t: 

Lemma 4.8. For general t, the solution xq to (|4.1(i|) has an expansion of the form 
1)4. 18(1 - 1)4. 19|) . The series in 1)4. 18(1 is absolutely and uniformly convergent for small 
e > 0, and when rearranged to 



(4.20) 



00 



1 s=0 



log log - 

e 



log- 

e 
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provides a complete asymptotic expansion for (xq/xq) — 1. The lead term in ()4.2U|) has 
coefficient doi = {t — l)^/r/. 

Proof. We give only a sketch of the proof. We can write the result of Lemma l4.3l in the 
form 

(4.21) w := -r/logxo = -a^^ - {t - l)cr"V - \ -v, 

with A := log{at-i/r]^~^), where v = o(l) and 

1 loglog(l/e) 
a := — : — — and r :— 



V 



log(l/e) ■ log(l/e) 

Substituting 1)4. 21(1 into the defining equation 1)4.16(1 for xq, we find 

t-1 t-1 
e = e'"Y,ak {-w/r^f = e"' r?"^ [a-^ + {t - l)a-\ + A + 

fc=0 A:=0 

or equivalently that w is a root of the function 

F(C) := - 1 - {[1 + (t - l)r + Act + ^C]*"^ - l} 

t-2 



A;=0 



The proof now proceeds as on pp. 347-349 in to obtain a series representation 
for f , and then we exponentiate to obtain ()4.18|) . We omit the details and will be content 
here to verify the asserted value of doi- Indeed, our proof sketch shows that rj doi equals 



2vrz 71(^1=^ 



-2 



where the integration is taken counterclockwise and G{Q is the coefficient of a^r^ in 
CF'(C)/F(C), namely, 

G(C) = {t- ine^ (e^ - l] 

Evaluating the integral completes the proof. □ 

Remark 4.9. In Section [S] we will see that, when viewed as an asymptotic expansion 
for xq, ()4.18|) is not accurate enough to use in obtaining a strong small-deviations result 
but is good enough (see Theorem 15. 3|) to yield an expansion for the logarithm of the 
small-ball probability. 

5 Small deviation estimates 

We now have at hand all of the ingredients we need to apply Sytaja's theorem and 
obtain p = p{e) := P{V'^ < e) up to a factor I + o{l), i.e., its log up to an additive o(l). 
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Until Section f5.3l we will assume d > 2. (The results for d = 1 are simpler to derive but 
require some modification.) 

To see how we are now positioned to obtain — logp up to additive o(l), note by 
Sytaja's theorem and Corollarv 14.41 that 

(5.1) - logp = h{x*) - ex* + i log [2^(1 - (i)exo] + o(l) 

where = l/(2mi +2) and xq is given explicitly by (|4.2|) . Here x* solves h'{x) = e, and 
Proposition 14 . 51 gives an asymptotic expansion allowing its computation to an arbitrarily 
large power of e once the solution xq to (|4.1j) is obtained (as discussed in Section . 
The expansion for x* can then be substituted into the expansion 1)3. 12() for /i(-). 

In this section we present for the reader's convenience some explicit small-deviations 
estimates in a few cases where the final result of the above program is reasonably clean. 
In Section [5. II we determine lead-order asymptotics, and more, for — logp. In Section l5^ 
we determine lead-order asymptotics for p itself. For completeness, in Section 15.31 we 
handle the one-dimensional case of m-integrated Brownian motion. 



5.1 Logarithmic small-deviations estimates {d > 2) 

Our first two results give lead-order asymptotics for — logp. 

Theorem 5.1. Let V be given as at (|1.4|) for an arbitrary {rrii, ... ,md) -integrated 
Brownian sheet with d > 2, and assume m := mi = ■ ■ ■ = nit < f^t+i < • • • < md- Then 
as £ ^ 0+ we have 2m+2 

2m + 1 rC(2m + 2)*-2 



-logP(y2 <e) = (l + o(l)) 



(2m + 1) 



t-i 



(t-l){2m + 2) 

2m + l / \\2m + \ / \\ 2m + l 

, log ■ 

el \ £, 



where 
(5.2) 



C(d,t,m) = C :-- 



sm 



\ 2m+2 j 



njr-=t+l ^rrij 



1 



2m+2 



(t- l)![(2m + 2)7r]'^-i 
Proof. Using Sytaja's Tauberian Theorem and Lemma 13.121 we find 

-logP(U2 <e) = (l + o(l)) [h{x*)-ex*] = [1 + o(l))2^-\l- i)ci,t-i {x*f{\ogx*)'~\ 

where ^ = l/(2m + 2) and ci^t-i is given by 1)3. 14() and so equals C. Now use Lemma 
and rearrange to obtain the desired result. □ 



Corollary 5.2. If mi = 

-logP {V^ < e) 

= (l+o(l))^ 



md = m, then as e ^ 0+ we have 



(2m + 2)-icsc 



2m+2 



(d-iy. [(2m + 1)it' 



d-l 



2m+2 
2m + l 



\\^2m + l 

£ 



log- 



(d-l)(2m+2) 
2m+l 



□ 
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Corollary 15.21 was previously established by Li . More recently, Theorem 15.11 
was obtained independently of us, and essentially simultaneously, by Karol' et. al [S] 
with a different form for the leading constant. Indeed, the trace sums that appear in 
Corollary 5.2 of jH] are simply generalized Dirichlet series values: recall our Remark l3.14l 

Next we show that it is possible to obtain arbitrarily high-order corrections to the 
estimate for —logp in Theorem 15.11 [Note, however, that the asymptotic scale clearly 
is not fine enough to obtain — logp up to additive o(l).] 

Theorem 5.3. The o(l) expression in Theorem, \5.1\ ho,s a complete asymptotic expan- 
sion of the form 



1 



log- 



1 



^"^Drsi log log ■ 

r=l s=0 ^ 

The lead-order coefficient is Du = {t — l)^(2m + 2)/(2m + 1). 

Proof. We give a sketch. Using Proposition 14.51 with j = 1, one finds that substitution 
of xq for X* in (|5.H) produces an error which is negligible relative to the scale involved 
in the statement of the theorem. Substitution of the expansion in Lemma 14.81 for xq 
leads to the desired result. The details are straightforward but quite tedious, and are 
omitted. □ 

Theorem 15.31 immediatelv gives us the following two examples. 

Example 5.4. When m = 0, Theorem 15.31 sharpens Csaki's classical result jH] for 
d-dimensional Brownian sheet (d > 2) to 



logP(y2 < e) 



(d-iy.TT' 



d-l 



-2 1 



log- 



2{d-l) 



l + 2(d- 1) 



log log ^ 
log^ 



+ 



1 



log- 



Example 5.5. For the (m, m)-integrated Brownian sheet, m > 0, we have 

2m + 2 



1 \ 2m + l 

a I 



I \ 2m+l 

log - + a „ , , 

e J 2m + 1 



1 1\2™+1 1 

-log - log log - 

e e £ 




where 
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5.2 Strong small-deviation estimates {d > 2) 

We focus on the case of distinct m^'s. In this case {t = 1) the statement of Theorem 15. II 
simphfies to -logP (F^ < = (i + o{l))E{e), where C is given by H5.2|l and E{e) is 
the expression 

2m + 2 1 
2m + l / C \2m + l / l\2m + l 

(5.3) E{e) = ^- 



2m + 2J \£ 

Theorem 5.6. Let V he given as at 1)1.4(1 for an {nii, ... ,171,1) -integrated Brownian 
sheet with d>2, and suppose that m := mi < • • • < m^. Then as e ^ 0+ we have 



P(^'<e) = (l + o(l)) 



m + 1 



-1/2 



exp{-£;(e) X [1 + S(e)]} 



for some finite linear combination S(e) of powers of e, wherein each power is a nonzero 
nonnegative integer combination of the numbers 

(2mi + 2)(2m^ + l) ^ 

(2mi + l)(2m^ + 2) ' " ' 

Proof. This fohows routinely from ()5.1|) and Example 14.61 recalling the notation r/,y = 
{2mi, + 1)/ (2mi, + 2) for = 1, . . . , d. We omit the details. □ 

Example 5.7. Even when d = 2, some vectors m require many terms in 5](e) in order 
to obtain the small-ball probability up to a factor 1 -|- o(l). We will be content here to 
illustrate Theorem 15.61 bv using it to state a first-order correction to Theorem 15.11 when 
d = 2 and mi < m2. The present example thus complements Example 15.51 
For the (mi, m2)-integrated Brownian sheet with mi < m2, we have 

1 2mi + 2 / 1 2mi+2- 



n /l\2"il+l / l\2rn-^+l' 2rn2 + 2 / / X \ 2mi +1 ' 2m2 +2 ' 

-l0gP(V^<£)=O,(-j +0,(-j +„^(- 

2mi+2 1 2mi+2 

where ai = (2mi -h l)[ci,o(l)]^^^, 02 = (2m2 + 2) [ci,o(l)]^^^'^^^C2,o(l), and 
ci^o(l) and C2,o(l) are given as in Lemma [3. 121 

Example 5.8. Suppose, as in Example 14.71 that m = mi = and that the m^'s are 
distinct, in which case (|5.2|) reduces to 

d 

C = {2n)-('-'^l[L*^^{-l/2) 
i=2 

and (|5.3|1 simplifies further to 

E(e) = — X -. 
^ ' 8 e 

Then, using the calculations in Example 14.71 it is easy to check that the only powers 
of e appearing in S(e) of Theorem 15.61 are powers 2r]i, — 1, for v = 2, . . . ,d, and 4r/2 — 2; 
we recall rj^ := (2mjy -|- l)/(2m,y -|- 2). If, in particular, d = 2, then only two terms are 
needed in S(e): see the next example. 
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Example 5.9. For (0, m2)-integrated Brownian sheet with m2 > 0, we have 

P(V^<e) = (l + o(l))^-^exp|... - + +.3(-) I 

1 2m2 

where di = [ci,o(l)]2, ^2 = (2m2 + 2)[ci,o(l)]^C2,o, and ds = [ci,o(l)]~^ [c2,o(l)]2, 
and ci^o(l) and C2,o(l) are given as in Lemma [3. 121 



5.3 Strong small deviation estimates when d = 1 



The difference between the cases d > 2 and d = 1 can be seen by comparing H3.16() 
and Lemma 13.21 and noting the constant term 2 log | detUm\ — (m + 1) log (2m + 2) in 
Lemma f3.2l Once this difference is noted, it is simple to obtain results for d = 1 like those 
in the preceding two subsections. In particular, the logarithmic small-deviations result 
in Corollarv 15.21 remains true verbatim when d = 1. Moreover, we have the following 
strong small-deviations analogue, studied previously in 0, 0, 0, ^Hj, and jH], of 
Theorem 15.61 

Theorem 5.10. Let V'^ = := £X^{t)dt, where Xm is m-times integrated Brown- 
ian motion defined in Remark M.ll Then as e ^ 0+ we have 



PiV^ < e) 



(l + o(l)) [|det?7„| •(2m + 2)("+i)/2 



vr 



m + 1 



E{e) 



'1/2 



exp{-S(e)}. 



Here, with C = csc(7r/(2m + 2)), the expression E{e) 
Vandermonde matrix defined at 1)3. 3() . 



IS 



given by (|5.3j) . and Um is the 
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